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基于 最 大 互相 关 UKF 的 传 感 网 目标 状态 和 
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摘 要 : 针对 空 天 地 一 体 化 传 感 网 络 中 传感器 观测 目标 时 观测 噪声 具有 重 尾 或 突变 性 质 的 问题 以 及 系统 偏差 对 目标 状 
态 估计 的 影响 , 提出 一 种 基于 最 大 互相 关 粒 无 迹 卡尔 受 滤波 (MCUKF) 的 目标 状态 和 系统 偏差 联合 估计 (ASMCUKF) 
算法 。MCUKF 算法 首先 通过 无 迹 变换 (UT) 获 得 预测 状态 估计 值 和 协 方差 矩阵 , 然后 使 用 基于 最 大 互相 关 羌 准则 (MCC) 
的 非 线性 回归 方法 重新 构建 观测 信息 ， 增 强 了 UKF 对 重 尾 噪声 的 鲁 棒 性 。ASMCUKF 算法 通过 目标 状态 向 量 扩 维 的 
方法 建立 状态 方程 和 带 有 系统 误差 的 非 线 性 观测 方程 ， 根 据 估 计 的 系统 偏差 进行 偏差 配 准 ， 改 善 了 系统 偏差 对 目标 状 
态 估计 的 影响 。 仿 真 结 果 表 明 ，ASMCUKEF 在 重 尾 非 高 斯 观测 噪声 的 环境 下 对 通信 目标 状态 和 系统 偏差 的 估计 效果 比 

传统 方法 更 好 。 
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Augmented target state estimation algorithm with systematic errors based on 
maximum correntropy unscented Kalman filter in sensor networks 
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Abstract: This letter propose a method that combines target state and system error based on the maximum correntropy criterion 
unscented Kalman filter( MCUKF), namely ASMCUKF, which can solve the problems that the observation noise of sensor is 
heavy-tailed or has some sudden change in Space-Air-Ground integrated sensor network. . Firstly, MCUKF uses unscented 
transformation (UT) to obtain the predicted state estimation and covariance matrix, and then a non-linear regression model based 
on MCC is used to reconstruct the observed information and strengthen the robustness of UKF in heavy-tailed noise. The 
ASMCUKF algorithm establishes the state equation and the nonlinear observation equation with systematic errors through the 
expansion of target state vector, and makes the error registration according to the estimated system errors to improve the 
influence of the System errors on the target state estimation. The simulation results show that the ASMCUKF has a better 
performance for the state estimation of the communication target than the traditional method in the environment of heavy-tailed 
non- Gaussian noise. 
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空 天 地 一 体 化 传 感 网 络 是 由 遥感 卫星 、 无 人 机 、 飞 艇 、 气 。 ”飞机 通信 、 车 载 网 、 工 业 控 制 、 环 境 监测 、 健 康 护理 等 领域 有 


球 、 地 面 传感器 等 传 感 资源 和 地 面 数据 处 理 中 心 组 成 的 新 型 网 ”广泛 的 应 用 前 景 上 3 。 
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录用 稿 


传 感 网 络 中 对 目标 相关 状态 估计 的 问题 ， 是 目前 研究 热点 
之 一 ， 主 要 通过 实时 的 优化 估计 算法 对 通信 目标 的 位 置信 息 进 
行 估计 预测 。 对 于 位 置 估计 常用 的 几 种 观测 方法 9 有 : 到 达 角 
(AOA) 、 到 达 时 间 (TOA) 、 到 达 时 间 差 (DTOA) 、 接 收 


的 信号 强度 (RSSI) ， 以 及 一 些 混 合 的 方法 。 因 为 目标 的 观测 
通常 是 非 线性 函数 ， 传 统 的 状态 估计 方法 主要 利用 卡尔 曼 相 关 


方法 对 通信 目标 状态 进行 估计 ， 并 假设 观测 噪声 为 高 斯 噪声 ， 
如 拓展 卡尔 曼 滤波 (EKF) 向 、 无 迹 卡尔 曼 滤 波 CUKF ) 器 ， 以 
及 一 些 改进 的 滤波 器 : 基于 最 大 似 然 估 计 的 卡尔 曼 滤波 四 、 基 
于 UKF 的 无 迹 粒 子 滤 波 CUPF)D。 而 在 传 感 网 的 实际 应 用 中 ， 
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关 烂 
六 OCD= 广 六 cr(a) (4) 
其 中 :6@= 坟 一 y，{(%,y))， 为 联合 分 布 函数 Fxy 的 M 个 采 


样 点 。 


2 MCUKEF 目标 状态 估计 算法 


MCUKF 算法 采用 基于 MCC 的 非 线 性 回归 方法 重新 构建 
见 测 信息 , 增强 了 UKF 对 重 尾 噪声 的 鲁 棒 性 , 文章 将 该 算法 应 
] 于 传 感 网 对 目标 的 状态 估计 中 ， 从 而 解决 了 传 感 网 中 重 尾 非 


~ 


Sag 


重 尾 非 高 斯 噪声 比较 常见 ， 而 EKF 和 UKF 都 是 基于 最 小 均 方 
差 准则 的 ， 对 重 尾 噪声 表现 出 强 敏 感性 ， 因 此 不 能 表现 出 良好 
的 滤波 性 能 外 。 为 了 解决 这 个 问题 ， 文 献 [9] 提 出 了 一 种 鲁 棒 的 
方法 一 一 互相 关 简 〈Correntropy) ， 其 不 仅 可 以 获得 二 阶 统计 
量 ， 还 可 以 获得 更 高 阶 的 统计 量 ， 基 于 它 的 准则 被 称 为 最 大 互 
相 闫 准则 (MCC) ， 并 且 已 在 一 些 领域 得 到 了 广泛 的 应 用 0 
1。 针 对 观测 噪声 呈现 重 尾 的 非 高 斯 状况 , 文献 [15] 提 出 了 一 种 
新 颖 的 MCUKEF 算法 ， 将 MCC 与 UKF 结合 提高 了 UKEF 的 鲁 
棒 性 。 
针对 单个 传感器 观测 的 不 确定 性 和 局 限 性 ， 通 常 利用 多 传 
感 器 进行 状态 估计 ， 但 其 没有 考虑 系统 偏差 带 来 的 影响 ， 容 易 
导致 估计 性 能 差 ， 甚 至 不 如 单个 传感器 估计 效果 明显 等 问题 ， 
文献 [16] 提 出 一 种 ASUKF 的 联合 估计 算法 , 对 系统 偏差 和 目标 
状态 进行 很 好 的 联合 估计 。 
本 文 将 MCUKF 算法 引入 到 传 感 网 中 ， 并 以 距离 和 角度 作 
为 传 感 节点 的 观测 量 ， 提 出 一 种 基于 MCUKEF 的 目标 状态 和 系 
统 偏差 联合 估计 算法 (ASMCUKF ) 。 在 考虑 系统 偏差 和 重 尾 
非 高 斯 噪声 的 情况 其 利用 MCC 对 非 高 斯 噪声 进行 鲁 
处 理 ， 并 通过 状态 向 量 扩 维 的 方法 来 估计 系统 偏差 ， 提 高 了 对 
目标 状态 估计 的 性 能 。 根据 核 宽度 取 值 不 同 , 将 ASMCUKF 算 
法 与 几 种 典型 的 联合 卡尔 曼 滤 波 算法 进行 性 能 对 比 ， 最 后 给 出 
仿真 结果 和 分 析 。 
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1 ”互相 关 灶 


定义 随机 变量 式 和 了 之 间 的 互相 关 焙 表示 为 
V(X,7)= Ex,(X -7 0 
其 中 : E 表示 期 望 算 子 ， Kx () 表示 核 函 数 ， o >0 为 核 宽度 ， 
通常 选用 的 核 函数 为 高 斯 核 函 数 : 


G- 力 ) (2) 


207 
假设 随机 变量 了 和 了 的 联合 分 布 函数 表示 为 ,x y) ， 则 
互相 关 焙 表示 为 


1 
K(X—y)=G,(x—y)= ew 


V(X,7)=|G,(x— WadFy (Gy) G) 


Hr 


对 传 感 网 络 实际 应 用 中 得 到 的 数据 量 有 限 ， 联 合 分 布 
Fxy 函数 未 知 的 情况 ， 可 通过 有 限 样本 平均 估计 量 来 估计 互相 


斯 噪声 影响 的 问题 。 
2.1 MCUKF 算法 描述 
定义 非 线性 离散 时 间 系 统 模 型 为 
X(E+D = fk, x(k)) + 9(k) (5) 
ykE+1) = h(xX(K+1D)) + rk+1) (6) 
其 中 : X(K+D 为 目标 的 状态 向 量 ， y(k +]) 为 传 感 网 络 的 观测 
向 量 ,f(.) 和 h(.) 分 别 表示 系统 和 观测 的 非 线性 函数 , 并 且 假 设 
它们 是 连续 可 微 的 ，9( 为 过 程 白 噪 声 , 其 协 方差 矩阵 为 Q(K) ， 
r(E+D 为 观测 白 噪 声 ， 其 协 方差 矩阵 为 RUK+D 。 
MCUKEF 滤波 器 有 具体 推导 如 下 05: 
首先 ， 根 据 UT 变换 获得 一 组 sigma 点 集 


也 


和 (KE 和 =XK| 有 
Zk]D= XE + (Yar DPEID) ,i=1~n (7 
ZEND=XEID (Vat DPEID) ,i=n+1~ 2n 


其 中 : (Var DPEID) 表示 和 矩阵 方 阵 的 第 


A=a?n+e)-n， Q 是 采样 点 的 分 布 函数 控制 参数 ， 一 般 为 
104<w<1，z 为 待 选 参数 ， 通 常设 置 为 3_m 。 
然后 根据 系统 方程 形成 2n +1 个 sigma 点 集 


i 列 ， 


XK+IR = (KY Kk|D),i=0~2n (8) 


则 系统 状态 量 的 一 步 预测 及 协 方差 矩阵 为 


KELL 各 = or +1 (9) 


i=0 


P(E+1|K) =7 {0 [Xk+1| 人 各 一 YOCE+I GD]x 
i=0 


名 (10) 
[XCK+1|K) -x K+ RT }+ QR) 
其 中 :采样 点 对 应 的 权 值 为 
wo =A/(n+A) 
@ =AMn+t+A)+(+oa +p) (11) 


: ; 1 
CO 3 0. 于 a 
2(n+1) 


B 为 非 负 的 权 系数 ,在 高 斯 情况 下 , 8=2 性 能 最 优 ,然后 根据 上 
一 步 预测 值 ,产生 新 的 sigma 点 外 


i=1~ 2n 


7 


录用 稿 


zx K+1| k= +1|R) 
YEH =XE+ID+ (Yat DPE+IAD) ,i=1~ n 
YEH = (Var DPE+ID) ,i=n+ 


i 


eA 人 /三 二 口 工 | 
maX /会 下 E 昌 日 干 | 
\ FA 


Ch 11ICVAAIVEC ] 门 | 
赵 季 红 ， 等 : 基于 最 大 互相 关 焙 UKF 的 传 感 网 目标 状态 和 系统 偏差 联合 估计 算法 


1 ~ 27 


(13) 


将 预测 的 新 sigma 点 集 代入 观测 方程 ， 得 到 相应 的 观测 


sigma 点 集 


Zi(k+1|K)=h Yt KR),i=1~ 2n 


(14) 


将 式 (5)(6) 等 式 结合 式 (9)(10) 重 新 得 到 以 下 的 非 线 怕 


X(k+1|k) X(k+1) n(X(k+1)) 
y(k+1) | PKE+LXKE+D) r(k+l) 


FE 模 型 : 
(15) 


对 等 式 (10) 对 应 的 协 方差 矩阵 PCk +1|k) 和 RUK+D 的 联 


合 矩 阵 y(k +1) 进行 Cholesky 分 解 
P(E+1|k) 0 

0 R(k+1) 
Ne 对 式 (15) 两 边 左 乘 T1(k +1) ， 得 到 如 下 等 式 : 
| D(k+1)= g(x(k+1))+e(k+1) 


ee | -Terr 


其 中 : 


D+D)=Tk+ oe | "| 


y(k+1) 


ser -Tt EY | 


h(x(k +1)) 


11(X(K 十 


= 
e(k+1)=T wa FED 


e(KE+D 中 第 i 行 的 元 素 e(E+D 为 
e(k+D)=d,(k+D)— g,(X(k+1)) 
定义 基于 互相 关 焙 的 开销 函数 为 


J(K+D = Yo +1)— g,(X(k+1))) 


其 中 : 4.() 、 g,() 分 别 为 向 量 D(.) 和 g() 的 第 i 个 元 素 
表示 向 量 D(.) 的 维 数 。 


最 优 估 计 ， 得 到 修改 后 的 协 方差 为 


) (10) 


(17) 


(18) 


(20) 


CD) 


， m+n 


基于 互相 关 灼 性 质 ， 当 函数 J() 最 大 时 可 以 得 到 状态 的 


DKE+D=T(E+DO (K+DT (K+ D=diag(Y.(k +D),Y,(k+D) 
(22) 
其 中 : 
CIKE+D=diag(G (e+D Ge (KE+D) (23) 
在 实际 中 真实 的 状态 是 未 知 的 ， 设 7()=0， 即 式 (15) 中 
XKE+D= 允 KE+I DO ， 则 有 
DK+D = P+1|R) (24) 
修改 后 的 观测 协 方差 为 
PK+D= 吃 (+D (25) 
则 系统 观测 的 先 验 均值 和 协 方差 为 
yk+IK) = VY wiz (er) (20) 
i=0 


2n 
A AVACERIL A 
= 1 吕 - KE+11D] 本 


[zs+1D-YE+HD] + Pr) 


互 协 方差 矩阵 和 卡尔 曼 增益 矩阵 为 


其 +LEH1 


>w[zdk AA 
i=0 


(28) 


K(k+D=P ,Pr 


XE 
KK+1JAK+L ?RHIYK4T 


(29) 


最 后 ， 计 算 状 态 更 新 和 协 方 差 更 新 
XKE+IKETD=XK+1|R) + KE+D YE+D -YE+1R) | 


G30) 
K(kK+1) B31) 


PUk+1E+D=PK+ID-KKE+DP 


HH 


3 ASMCUKF 目标 状态 估计 算法 


ASMCUKEF 算法 将 系统 偏差 作为 状态 向 量 中 的 分 量 ， 利 用 
MCUKEF 滤波 器 实现 目标 状态 和 系统 偏差 的 联合 估计 ， 在 重 尾 
非 高 斯 观测 噪声 情况 下 ， 解 决 了 多 传感器 系统 偏差 对 通信 目标 
状态 估计 性 能 影响 的 问 
3.1 系统 描述 

空 天 地 一 体 化 传 感 网 络 中 ， 部 署 工 个 传感器 对 目标 状态 进 
行 估计 ， 传 感 器 的 坐标 为 (x,, y。,, zo,) 。 假 设 移动 通信 目标 第 k 
时 刻 的 位 置 和 速度 为 (x(), y(R),z(K)) 、 (Gy ， 传 感 
器 的 观测 量 为 (q,(1),8(k),g(k)) ; 工 个 传感器 的 距离 偏差 、 水 平 
偏向 角 偏 差 、 仿 仰角 偏差 为 x (Kk)=(Ad,,A9,Ap) ， 其 中 
i=7,2,.…, 上 。 则 定义 时 刻 目 标 状态 向 量 为 

X (Kk) = x(k), Xk), yk), YR), (Kk), kK) G32) 

假设 目标 做 匀速 直线 运动 ， 系 统 的 状态 将 由 目标 状态 和 系 

统 偏 差 的 联合 组 成 ， 则 扩 维 后 的 状态 方程 为 


离 “ 
o 


Xx(k+1)=F. Xx(k)+G.wk) (33) 
其 中 : 
X(K)=[X, (Kk) ,Xx (Kk) ] G4) 
W(k) =[q,,43,93] 
F=diag(F.,)),G= diag(G,,0) (35) 
[1 T0000 I 0 0| 
0 10 0 0 0 T 0 0 
0 0 17 0 0 2 
请 = ,G, = 0 0 (36) 
0 0 0 10 0 0 了 0 
0000 17 0 0 三 
[0 00001| 0 0 7 
其 中 了 为 单位 阵 。 
系统 的 观测 方程 定义 为 
Z(K) = h(X(K)) + Ae(K) +Y(K) G37) 
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其 中 : 
h(X(K))=[d;(k), 0.(K), p(k)]= 
(XR) — X00) + (YR) — yo) + zk) —z0)" 
PR 
X(k) 一 Xi 
arctan US Bk 
| (xX(K)— X00) + VK) yo0) + (zk) 一 zo7 
其 中 Ae(k) = (Ad,,A0,Agp,) 表示 第 i 传感器 大 时刻 观 测 时 对 应 

的 系统 偏差 。 


3.2 ASMCUKF 算法 描述 

该 算法 是 以 MCUKF 算法 为 基础 ， 根 据 式 (33) 状 态 方程 和 
式 (37) 观 测 方程 来 配 准 系统 偏差 ， 并 改善 重 尾 非 高 斯 噪声 对 目 
标 状态 估计 的 影响 。 具 体 算法 步 又 如 下 : 

a) 选 择 合适 的 核 宽度 o， 设 置 初 始 的 状态 值 X(0|0) 和 协 方 
差 矩 阵 P(0|0) ， 时 刻 上 =1; 


I 


b) 根 据 式 (7) 计 算 sigma 点 ， 并 通过 扩 维 后 的 状态 转移 矩阵 
五 计算 出 sigma 点 集 的 一 步 预测 ; 
XI K+IK)=F:. Yk|R,i=0~2n (38) 
co) 通 过 式 (9)(10) 得 到 扩 维 后 状态 向 量 的 一 步 预测 
KE+I 有 及 协 方差 矩阵 PE+1|AD ; 
PTD=Y oT XK xO] 
气 (39) 
[RE+H 有 -YOU+ID] +G.Q0D:Gr 


d) 利 用 式 (13)(14) 产 生 新 的 sigma 点 集 和 对 应 的 带 有 系统 仿 
差 的 观测 值 ; 
e 通 过 式 (15)~(25) 得 到 修改 后 的 观测 协 方差 矩阵 请 KK+D， 
j 式 (18) 计 算 观 测 的 先 验 均值 六 K+lk) 和 协 方差 已， ; 

全 用 式 (30)(31) 得 到 对 系统 偏差 进行 配 准 的 状态 估计 和 协 


是 当 大 等 于 预先 设 定时 间 步 长 时 , 算法 结 
步骤 b)。 


4 ”仿真 分 析 


4.1 仿真 环境 
本 仿真 采用 3.1 节 的 系统 来 对 
RMSD 以 及 TARMSD 为 基准 进行 性 能 比较 ， 


, 否则 Kk=kK+7， 


涡 


标的 状态 进行 估计 ， 并 以 


1 M 
RMSD(K) = (3 E(x(K) -EK|R)) + 


m=1 


1 (40) 
OE) -HEL RY + (zk) — 2k] A)) I 
大 一] 天 
TARMSD = 二 六 RMSD( (41) 


其 中 :M=100 表示 蒙特 卡 罗 实 验 的 总 次 数 ，K=500 为 每 次 蒙特 
卡 罗 实 验 中 的 整个 时 间 步 长 。 
目标 在 传 感 区 域内 做 匀速 直线 运动 [17]， 运 动 轨迹 为 


Chi 
红 ， 等 : 基于 最 大 互相 关 焙 UKF 的 传 感 网 Bl 村 计算 法 


体 表 述 形式 如 下 : 


maX 上 iv 合 作 期 


日 二 || 
| 


x(k) = 6000 + SOk 
y(k) = 50000—10k 
z(k) =1000+8k 


在 传 感 区 域内 部 署 四 个 传感器 ， 其 坐标 和 对 应 的 系统 偏差 
分 别 为 (距离 单位 为 m, 角 度 单 位 为 rad) 


A(0,40000,0), B(40000,40000,0) 
C(0,60000,0),D(40000,60000,0) 

Ad = 2000,Ad, =1500,Ad, =1500, Ad, = 2000 

Ag = 0.0087,A@, = 0.0025,Ag = 0.0070, A@, = 0.0078 
Am = 0.0087,Ap, = 0.0025,Am, = 0.0070,Am, = 0.0078 


观测 精度 均 为 
0O0=05=03=004=1007， 
oo=0o=0o=oo=0.5 ， 
O01=0,,=0%3=0%4=0.5 


= 
设 实验 中 目标 状态 向 量 如 式 (34) 所 示 ， 初 始 真实 状态 为 (位 
置 单位 m, 速 度 单位 m/s) 
X(0|0) = (6000,50,50000,-10,1000,8 ,0,0,0)” 
初始 状态 估计 值 和 协 方差 矩阵 为 
X(0|0) = (6250,50,50500,-10,1250,8,500,0.001,0.001D)7 
P(0|0)= diag([100,2,100,2,100,2,10°,6x10°,6x10°]) 
过 程 噪声 为 
~ N(0,10m’) 0 
观测 的 非 高 斯 噪声 为 
] 一 0.8N(0,(0270)2)+0.2N(0,(10370D2) 


~ N(0,10m’) gq; ~ N(0,10m?) 


ys ~ 0.8N(0,((0.5z /180)raqd)’)+0.2N(0,((S5z /180)rad)’) 

7, ~ 0.8N(0,((0.5z /180)rad)’)+0.2N(0,((Sz /180)rad)’) 
4.2 ”仿真 结果 

系统 观测 目标 时 ， 在 受到 重 尾 非 高 斯 噪声 的 影响 下 ， 对 本 
文 算 法 ASMCUK 与 ASUKF00 和 ASEKF07 性 能 进行 了 比较 。 
图 1 为 通信 目标 在 传 感 网 中 移动 时 各 种 滤波 器 的 估计 轨迹 。 
在 处 理 基 于 非 线 性 方程 的 观测 量 时 , ASUKF 和 ASMCUK 显然 
要 比 ASEKF 估计 更 加 精确 ， 在 处 理 非 高 斯 的 观测 噪声 时 ， 


< 


ASMCUK 算法 要 比 ASUKF 更 好 的 鲁 棒 性 。 显而易见 ， 由 
于 初始 的 估计 与 真实 的 状态 具有 一 定 的 偏差 ， 所 以 估计 逐渐 通 
近 真 实 值 。 
6000 
© sensor 


real path 


Da 


1 各 种 滤波 器 对 目标 的 轨迹 估计 


十 | 


49 工 i| 
门 
3 起 所 半身 法 


图 2 


和 述 了 传 感 网 络 中 通信 目标 在 不 同 滤波 器 下 的 RMSD。 
随 着 迭代 次 数 的 增加 ， 各 种 滤波 器 的 RMSD 呈 减 小 趋势 , 但 1 
于 观测 受到 重 尾 的 非 高 斯 噪声 影响 ,所 以 RMSD 曲线 出 现 明显 


ASUKF 进行 相 比 ， 新 算法 表现 出 了 对 重 尾 非 高 斯 噪声 的 强 鲁 
性 。 仿 真 结果 表明 ， 在 重 尾 的 非 高 斯 噪声 的 影响 下 ， 较 大 的 
核 宽度 (o =20) 时 ，ASMCUKEF 将 表现 出 与 ASUKF 类 似 的 性 


误 


波动 。 在 处 理 基于 非 线 性 方程 的 观测 量 时 ， 由 于 ASUKF 类 型 ” 能 ， 较 小 的 核 宽度 (oa =2) 时 ， 可 以 实现 比 ASEKF 和 ASUKF 
滤波 器 利用 无 迹 变换 (UT) 在 估计 点 附近 确定 采样 点 ， 用 这 些 。 更 好 的 性 能 ;并 且 随 着 传感器 的 增多 , 估计 性 能 也 能 得 到 有 效 的 
采样 点 表示 的 高 斯 密度 去 近似 表示 状态 的 概率 密度 函数 ， 因 此 ”改善 。 

表现 出 比 ASEKF 滤波 器 更 好 的 性 能 。 在 处 理 非 高 斯 的 观测 噪 

声 时 ， 由 于 互相 关 灼 强 鲁 棒 性 ，ASMCUKEF 在 所 有 的 滤波 器 和 二 


一 6 一 ASMCUKF(o=2) | ] 
一 全 一 ASMCUKF(c=3) 
ASMCUKF(a=5) 
一 全 一 ASMCUKF(c=10) | J 
一 + 一 ASMCUKF(c=20) 


中 RMSD 最 小 , 如 表 1 所 示 , 使 用 适中 的 核 宽度 时 , ASMCUKF S00 
明显 优 于 ASUKF。 特 别 o=2 时 ，ASMCUKEF 表现 出 最 小 的 
TARMSD， 核 宽度 非常 大 时 ，ASMCUKEF 的 性 能 与 ASUKF 基 


本 相同 。 


600 性 300 上 


a ASEKF 
"| | 本 : ; : : 
ASMCUKF(o=5) 传感器 个 数 /个 
天 由 ASMCUKF(o=10) | ee 
ASMCUKF(o=20) | 图 3 多 传感器 时 目标 估计 误差 
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